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OPTIMALITY OF REFRACTION STRATEGIES FOR SPECTRALLY NEGATIVE LEVY 

PROCESSES 


DANIEL HERNANDEZ-HERNANDEZ*, JOSE-LUIS PEREZ*, AND KAZUTOSHI YAMAZAKit 

Abstract. We revisit a stochastic control problem of optimally modifying the underlying spectrally nega¬ 
tive Levy process. A strategy must be absolutely continuous with respect to the Lebesgue measure, and the 
objective is to minimize the total costs of the running and controlling costs. Under the assumption that the 
running cost function is convex, we show the optimality of a refraction strategy. We also obtain the conver¬ 
gence of the optimal refraction strategies and the value functions, as the control set is enlarged, to those in the 
relaxed case without the absolutely continuous assumption. Numerical results are further given to confirm 
these analytical results. 
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1. Introduction 

We consider a class of stochastic control problems under the restriction that strategies must be absolutely 
continuous with respect to the Lebesgue measure. Given a spectrally negative Levy process, the controller 
aims to optimally decrease the process so as to minimize the sum of the expected running and controlling 
costs. The running cost is modeled as a convex function h of the controlled process that is accumulated 
over time. The controlling cost is proportional to the size of control. Our objective is to show the optimality 
of the refraction strategy so that the controlled process becomes a refracted Levy process of Kyprianou 
and Loeffen [15], with a suitable choice of the refraction boundary level. 

In the last decade, there have been great developments in the theory of Levy processes and its applica¬ 
tions in stochastic control. A representative example is its contributions in de Finetti’s optimal dividend 
problem, where the classical compound Poisson model has been generalized to a spectrally negative Levy 
model. The expected net present value (NPV) of total dividends under the reflection strategy, so that the 
controlled process becomes a reflected Levy process, can be concisely written in terms of the scale func¬ 
tion (see [2]). While the optimality may fail depending on the choice of the Levy measure, Loeffen [17] 
obtains a sufficient condition for optimality and in particular shows that it holds when the Levy measure 
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has a completely monotone density. Other stochastic control problems that have been explicitly solved via 
the scale function include the dual model of the optimal dividend problem [5], inventory control problems 
as in [4, 21], and games between two players [11]. 

While the reflection strategy is commonly shown to be optimal in these papers, implementing it in real 
life is often not feasible nor realistic. Most likely in every scenario, the ability of controlling is limited and 
hence is more reasonable to restrict the amount of modification one can make to the underlying process. 

In classical settings of the optimal dividend problem, reflection strategies that are shown to be optimal 
always lead to ruin in finite time; these results have been in dispute and there has been a need for new 
models so as to avoid these undesirable conclusions. In currency rate control (see, e.g., [12] and [18]), 
where a central bank controls the currency rate so as to prevent it from going too high or too low, reflection 
strategies, that are shown to be optimal under many models, are not implementable in real life. In 2011, 
Swiss central bank introduced the peg so as to stabilize their currency in response to its escalating value 
against other currencies: they set a ceiling on its rate against Euro, and literally implemented a reflection 
strategy. However, this only lasted until the beginning of 2015 when they had to admit the enormous cost 
of doing it and scraped the ceiling. 

One way to model these limitations on the control set is to add an extra restriction on the set of admissi¬ 
ble strategies to be absolutely continuous with respect to the Lebesgue measure. Kyprianou et al. [16], in 
particular, consider the optimal dividend problem under this restriction. They show under the completely 
monotone assumption on the Levy measure that a refraction strategy is optimal: Roughly speaking, the re¬ 
sulting controlled process becomes a refracted Levy process that progresses like the original Levy process 
below a chosen threshold while it does like a drift-changed process above it. 

In this paper, we revisit the stochastic control problem that minimizes the convex running costs. The 
setting is similar to the existing literature on singular/impulse control problems as in [6], [7] and [21]. The 
convex assumption is typically needed and is required in, e.g., [3, 9, 20]. Differently from these papers, we 
consider a version with an additional condition that the strategy is absolutely continuous. More precisely, 
a strategy must be of the form Lt = J^isds, t > 0 , with i restricted to take values in [0, 5] uniformly 
in time. As in [16], our objective is to show the optimality of a refraction strategy. While a reflected 
Levy process can be dealt relatively easily because it moves just like the original Levy process except at 
the times running maximum/minimum is updated, a refracted Levy process moves in two different ways 
depending on whether it is above or below the given threshold. The computation thus becomes more 
intricate. 

In order to tackle our problem, we take the following steps. 


(1) By directly using the resolvent measure for the refracted Levy process as in [15], we first write the 
expected NPV of total costs using the scale function under each refraction strategy. 

(2) The candidate threshold b* E [— cxo, cxo] is then chosen by observing the identity (3.20) below. We 
show that such b* is given as a root of a monotone function / defined in (3.23) below, or otherwise 
it is either — cxo or cx). 
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(3) We then verify the optimality of sueh strategy. Toward this end, we derive a suffieient eondition 
for optimality (verifieation lemma), and then show that it is equivalent to certain conditions on the 
smoothness and slope of the value function. We observe that the level b* is such that the candidate 
value function becomes smooth and convex; using these facts, we complete the proof of optimality. 

It is important to highlight the fact that this set of results are obtained without using directly the HJB 
equation associated with this control problem and, instead, it is used in an indirect way to derive sufficient 
conditions for optimality, based only on the derivative of the value function; see Lemmas 4.1 and 4.3. 

In addition to solving this problem, we analyze the behavior as the upper bound 5, of the control set, 
becomes large. Being consistent with our intuition, we show that the optimal refraction strategy of this 
problem converges to the reflection strategy that has been shown to be optimal in the original problem with¬ 
out the absolutely continuous condition. In particular, we show the convergence of the optimal threshold 
b* and value function to those obtained in [21]. In order to discuss the practical side of implementing 
the obtained optimal strategy, we give a series of computational experiments using the phase-type Levy 
process and quadratic cost function. We confirm the optimality as well as the convergence as 5 t oo to the 
reflection strategy. 

The rest of the paper is organized as follows. The problem is defined in Section 2. In Section 3, we first 
review the refracted spectrally negative Levy process and the scale function. We then compute the NPV 
of total costs under a refraction strategy and choose a candidate refraction boundary level b*. Section 4 
verifies the optimality of the selected refraction strategy. Section 5 studies the convergence as 5 t oo of 
the optimal refraction strategy to the optimal reflection strategy. We conclude the paper with numerical 
results in Section 6 . 


2. Mathematical Model 

Let (fl, P) be a probability space hosting a spectrally negative Levy process X = {Xp f > 0} whose 
Laplace exponent is given by 

2 r 

(2.1) i/;(6') := logE[e®^i] = 76*-f -f / - I - 6> > 0, 

J (—CXD, 0 ) 

where v is a Levy measure with support in (— 00 , 0) and satisfying the integrability condition q)(1 A 
z^)iy(dz) < 00 . It has paths of bounded variation if and only if a = 0 and | 2 ;| ly(dz) < 00 ; in this 

case, we write ( 2 . 1 ) as 

^(^e) = de+[ (e®" - i)z/(dz), e>o, 

J (—00,0) 

withy ;= 7 — 2 ; z/(dz). We exclude the case in which X is the negative of a subordinator(i.e., X has 
monotonically decreasing paths a.s.). This assumption implies that 7 > 0 when X is of bounded variation. 
Let Pa; be the conditional probability under which Xq = x (also let P = Pq), and let F := {Xp f > 0} be 
the filtration generated by X. 

Fix /3 e M, (5 > 0 and a measurable function h : M ^ M that is specified in Assumption 2.2 below. 
Define II^ as the set of absolutely continuous strategies I given by adapted processes L* = f* 4ds, t > 0, 
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with i restricted to take values in [0, 5] uniformly in time. For g > 0 fixed, the objective is to consider the 
NPV of the expected total costs 


( 2 . 2 ) 

where 


Vi[x) := Ej 




./o 


\J[ := Xt -U, t> 0 , 

and compute the (optimal) value function 

(2.3) v(x) := inf vAx), x G M, 

as well as the optimal strategy that attains it, if such a strategy exists. 

We make the following standing assumptions on the Levy process X and the running cost function h. 


Assumption 2.1. (1) For the case X is of bounded variation, we assume that 7 — 5 > 0. 

(2) We assume that there exists 9 > 0 such that exp(0|z|)z/(dz) < 00 . 

Assumption 2.2. We assume that h is convex and has at most polynomial growth in the tail. That is to 
say, there exist ki, /C2 > 0 and N eN such that h{x) < ki + k 2 \x\^ for all x G M such that |x| > m. 


Define the drift-changed Levy process 
(2.4) Yt := Xi -5t, t> 0, 


which is the resulting controlled process if (. is uniformly set to be the maximal value 6. Assumption 2.1 
(1) guarantees that Y is not the negative of a subordinator, and hence is again a spectrally negative Levy 
process. This is commonly assumed (see, e.g., [16]) so that the refracted Levy process given below is 
well-defined. 

By Assumption 2.1 (2), the domain of the Laplace exponent (2.1) can be extended to (—0, cxo) and, by 
its continuity, we can choose sufficiently small 6 > D such that f{—9) < q. With this choice of 6, 

poo poo 

(2.5) E[e-®^'=<^] = q e-''*E[e-®^‘]df = q < 00 , 

Jo Jo 

where Cq is an independent exponential random variable with parameter q. 

Assumptions 2.1 (2) and 2.2 guarantee that vi as in (2.2) is well-defined and finite for all i G II^. To see 
this, by the convexity of h and because Yt < < Xt, 

(2.6) h{Ui) <h{Xt)y h{Yt), 


while, with h ;= inf^-gK h{x) G [—cxo, cxo), we have 


h if h > — 00 , 

h{Xt) A h{Yt) otherwise. 


(2.7) 


h{Uf) > 
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Hence, for any random time T (independent of Cq), 


E.t 


( 2 . 8 ) 


sup \h{Ul)\dt 


< E, 


-'?‘(|/i(Xt)| + |M>"t)l + IZf|l{/i>-oo})df 


= q ^E^ 


Ue.>T}mXj\ + \h{Yj\ + \h\l 

{/),>—CX)}) 


< oo, 


where the finiteness holds by (2.5) and Assumption 2.2. The finiteness (2.8) will be particularly important 
in the verification step (Lemma 4.1 below). As is commonly used in verification arguments, localization 
is first needed; (2.8) lets us interchange the limits over integrals. 


Remark 2.1. We can also consider a version of this problem where a linear drift is added to the increments 
of X (as opposed to being subtracted): one wants to minimize, for some /3 G M, the NPV 


ve{x):=E, / e-^\h{Xt + L^) + pit)dt 


'-JO 


However, it is easy to verify that this problem is equivalent to the problem described above. 

To see this, we use Y as in (2.4) and set Lt := 5t — Lt (which is admissible with Lt = J^isds where 
£s '.= S — £s E [0,S] a.s.). Then, we can write 


ve(x) = 




ih{Y, 


^£t)dt 



q 


Hence it is equivalent to solving our problem for fd := —fd. 


3. Refraction Strategies 

One of the objectives of this paper is to show the optimality of the refraction strategies, say t G II^, 
under which the controlled process becomes the refracted Levy process = {U^]t > 0} of [15], with 
a suitable choice of the refraction boundary 6 G M. By [15, Theorem 1 and Remark 3], this is a strong 
Markov process given by the unique strong solution to the SDE 


dU^ = dXt-dl^jjb^^^dt, f > 0. 


Namely, progresses like X below the boundary b while it does like Y above b. Let us write the 
corresponding NPV of the total costs associated with £^ by 


(3.1) 


Vb(,x) 





X eR. 


3.1. Scale functions. As it has been studied in [15], the expected NPV (3.1) can be expressed in terms of 
the scale functions of the two spectrally negative Levy processes X and Y. Following the same notations 
as in [15], we use and for the scale functions of X and Y, respectively. Namely, these are the 
mappings from M to [0, cxo) that take value zero on the negative half-line, while on the positive half-line 
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they are strietly inereasing functions that are defined by their Laplace transforms: 

1 


—Ox 


(3.2) 


where 


W^^\x)(\x = 


.-exw(9)(2:)dx = 


'tp{e) - q 
1 


'ijj{9) — 69 — q 


, 9 > <f)(g), 

, 9>^{q), 


(3.3) $(g) := sup{A > 0 : -^(A) = q} and (p{q) := sup{A > 0 : -^(A) — 6\ = q}. 

By the strict convexity of we derive the strict inequality (^(g) > $(g) > 0. 

Scale functions appear in the vast majority of known identities concerning boundary crossing problems 
and related path decompositions. In order to illustrate the importance of the scale functions we provide 
the following result, the so-called two-sided exit problem (see for instance Theorem 8.1 in [14]). Define 

r/ = inf{t >{)■. Xt > a] and Tq = inf{f > 0 : < 0 }. 

Then for all g > 0 and x < a, (x)(a). For other applications of the scale 

function, see, among others, [13] and [14]. 

Fix A > 0 and define ^a(-) as the Laplace exponent of X under with the change of measure 

= exp{XXt —'ijj{X)t), t > 0] 

see [14, Section 8 ]. Suppose that is the scale function associated with X under P^ (or equivalently 
with '0 a('))- Then, by Lemma 8.4 of [14], x G M. In particular, we define 

(3.4) IF$(,)(a;) := kF®)(x) = a; G M, 

which is known to be an increasing function and, as in Lemma 3.3 of [13], 

(3.5) IF#(q)(a;) ^'0^(T’(g))~^ < cxo as a; —)■ cxo. 

Below, we summarize the properties of the scale function that will be necessary in deriving our results. 


dP^ 

liP 


Remark 3.1. (1) If X is of unbounded variation or the Levy measure is atomless, it is known that 

is C'^(M\{0}); see, e.g., [8, Theorem 3]. 

(2) Regarding the asymptotic behavior near zero, as in Lemmas 3.1 and 3.2 o/[13], 


I ^ if X is of unbounded variation, 
' 1 if X is of bounded variation. 


(3.6) 


2 


if a > 0 , 

pj/('?)'( 0 -)-) ;= lim IF^'^^'(x) = cxo if a = 0 and z/(—cx),0) = cxo, 

q+iy(-oo,o) if (j = D and i>{—oo,D) < oo. 


(3) Asin(8.22) and Lemma 8.2 of [14], < W^‘^^'{x+)/W^'^'>{x)fory > x > 0. In 

all cases, W^^^fx—) > {x4-) for all x > 0 . 
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We also define, for all x > 0, 

(3.7) := > 0. 

Let X_f := info<i'<t X^i, f > 0, be the running infimum proeess. By Corollary 2.2 of [13], for Borel subsets 
in the nonnegative half line, 

(3.8) P{-^eq e dx} = “ gW"('')(x)dx = ^[0^'^^(a:)dx + iy('')(0)(5o(dx)], 

where fL*^'?)(dx) is the measure sueh that PL('^)(x) = PL('^)(d^) (see [14, (8.20)]) and (5o is the Dirae 

measure at zero. 


Remark 3.2. (1) For the case of spectrally negative Levy process X, 



q 6 — $(g) 

$(g) V'( 0 ) - 



and Xeq is exponentially distributed with parameter $(g) where Xt := supo<t'<t Xt', t >0, is the 
running supremum process; see Section 8 o/[14]. 

(2) It is easy to see that 


0('')(M)dM = ^ 0} = ^ - iy(^)(o), 


(3.9) 


uQ^^\u)du= ^E[- 






$(g) 


E 


7 


ip(q)X 


p{q) - ^{q) 

8p{q) 


-PL(«)( 0 ). 


For any —M < x, by (3.7) and (3.8), 


/ —M p—M 

|/i(i/)|e-^('')W4(,)(x - y)dy = / |/i(i/)|0(^)(x - y)dy 

■oo J—oo 

= i " |/!(y)|P{2C„ + X e dy} = L«)e-»W-E,||ft(2i:^J|l,^^<_„,]. 

7 J —oo 7 

Now, reeall Assumption 2.1 (2) and fix any 0 < 6* < i9. By the independenee of X^^ — X„ ^ and X„^, 
E[g-^'Veq] _ whoro in partieular Xeq — ~ Xeq (by duality) is exponentially 

distributed with parameter $(g) as in Remark 3.2 (1). Henee, we see that 


(3.11) 


E[e®'^q'] < oo. 


By (3.11) and Assumption 2.2, (3.10) is finite. 
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3.2. The expected NPV (3.1) in terms of scale functions. Fix 6, x e M. By Theorem 6 (iv) of [15], the 
resolvent measure 


Rb{x,B) := e 5 } = E, 

admits a density 

(3.12) 
given by 

(3.13) 




, Be 


Rb{x,dy) = {rl^\x,y) + r^^\x,y)l{^>b})dy, yeR, 


rb\x,y) ■ = 


rb\x,y) ■ = 

where we define 

(3.14) M(x; &) := ((^(g) — <F(g))e“'*’^'^^^ y — 2;)dz, x > b. 

Henee, (3.1) ean be written (see Remark 3.3 below) as 

(3.15) Vb{x) = 'yi^^(x) + v^^\x)l{^yb}, 


tp(q)-^(q) ^^(q)(x-b)-ip(q)(y-b) V > h 

S^{q) ^ — ’ 

g$(g)(x-b) y(gKj&(g) gy(g)6 Q-<p(q)zy^(qy^^ _ - W^i^{x - y) y < h, 

e-‘^('^)(*'-^)M(x; b) — W^'^^(x — y) y > b, 

(5M(x;6)e^('')^/“e-^('?)W(«)'(z-|/)d^-5 W^i\x - z)W^'J^'{z - y)dz y < b, 


where 

(3.16) 


(3.17) 




:= e 


h{y + b) 


mq) 

^{q){x-h) y^i.q) - ^{q) 

$(g) 


h{y + + 


IL 

^(q) 


- y)dz - W^'^\x -b-y) 


dl/, 


vf\x) := ^ {h{y + b)+ /3(5)|e-‘^('')^M(x; b) - W^'^\x - b - y)^dy 

/ O /* oo /*x 

h{y + b)\^M{x;b) / e-^^'^^^W^'^^'{z-y)dz- W^‘^\x - z)W^‘^^’{z - b - y)dz^dy. 


Remark 3.3. In view of (3.15), both and are finite and hence the decomposition is well-defined. 

To see this notice that, by the estimation (2.8), v^^\x) is finite for x < b. Moreover, in view of the form of 
in (3.13), v\^\x) for X > b is also finite. Indeed, for any x > b and y < b, 

rl^\x,y) = e'*’d)(^-'>)rW(5^^) _ W^‘^\x - y) + e^d)(^-b)p^(9)(5 _ yy 


By the decomposition (95) of [\3], —W^’^\x — y)-\-e^^‘^">^^ ^">W^’^\b — y) = — ^^u{b — y) 

(where is the q-resolvent density of —X as in Theorem 2.4 (iv) o/[13]), with which h is integrable by 
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(2.8). Finally, using again (2.8), we conclude thatvj^\x) is finite as well, for x > b. In addition. 


(3.18) 


\h{y)\\r^^\x,y)\dy < oo. 


Since \h{y)\{rl^'* + r 2 ^^)(x, y)dy < oo and > 0, we must also have 


(3.19) 


\h{y)\\rf\x,y)\dy < oo. 


3.3. First order condition. We shall first obtain our candidate refraction boundary h*. In view of (3.16) 
and (3.17), (once it is confirmed that the derivatives d/dh and d/dx can be interchanged over the integrals) 
the next identity holds: 

d 


(3.20) 

where 


dh 


Lb(x) = Mb(x), xf-h. 


Uh{x) := 


-v[{x), x^h. 


'-Jo 


Since the first-order condition dvb{x)/db\b=b* = 0 is a necessary condition for the optimality of the 
refraction strategy , we shall pursue b* such that Ub*{x) vanishes; consequently, (if such b* exists) 


viAx 


= Ej. [ e . This identity will be important later in the verification of optimality. 


Let us define, for fc = 1,2, 


(3.21) := I h'{y)rl''’{x,y)dy-vl^’^x), x ^ b, 

J —OO 

so that 

(3.22) Ub{x) = + uf\x)l{:c>b}, X f^h. 

The next lemma shows that the right hand side of (3.21) can be written succinctly using the function: 

(3.23) 


Sk) 


(k)!, 


m ■■= 


Fig) - ^{q) 

Fig) 

Fig) - ^(g) 

Fig) 

5 ^ 


h\y + h)e-^^‘i^ydy + 5\ / h\y + h) / - y)dzdy - f 


h\y + h)e-^^'i'^ydy 


Fig) Lj- 


h'iy + 


<h(g)) 




^(g) 

Fig) 


mg) 


where the second equality holds because, by integration by parts. 


(3.24) 


-^^^m^^^'{z-y)dz + W^^\-y) = Fig) - y)dz, y<d. 


The proof of this result is technical and long (partly because we need to confirm that the derivative d/dx 
can go into the integral over the unbounded set), and therefore we shall defer the proof to Appendix A. 
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Lemma 3.1. We have 

(3.25) - h(h)W^‘^\x - b), x ^ b, 

(3.26) u^^\x) = {M{x]b)-W^‘^\x-b))I{b) + h{b)W^'^\x-b), x > b. 

In view of (3.22), this lemma direetly implies the following. 


Proposition 3.1. For all x, 6 G M such that x ^ b, 


(3.27) 


Ub{x) 


' FiQ) - ^(,)(x-b) 

[ sm 


+ (M(x; 6 ) - W('?)(a; - b))]l{b). 


3.4. Candidate value function. In view of our diseussion in the previous subseetion and Proposition 3.1, 
the ehoiee of our eandidate threshold level b* is elear. In (3.27), the braeket term on the right hand side is 
uniformly positive for all x G M. Indeed, for x > 6, it equals the resolvent density rl^\x, b) + rf\x^ b) 
at b in view of (3.13). Henee, we shall pursue b* sueh that I{b*) vanishes. However, the existenee and 
uniqueness of sueh b* is not guaranteed, and henee we shall define b* earefully here. 

First, by (3.7), the first equality of (3.23) and the eonvexity of h, the funetion I is nondecreasing. Henee 
we ean define the limits /(oo) := lim^^oo-^(^) ^ (— 00 , 00 ] and/(— 00 ) := \imi,i-oo Ii^) ^ [— 00 , 00 ). 
We set our eandidate optimal threshold level b* to be any root of /(&) = 0 if /(— 00 ) < 0 < /(oo). If 
/(oo) < 0 , we let b* = 00 ; if /(— 00 ) > 0 , we let b* = — 00 . 

In faet, if h is affine, then / is a eonstant. Conversely, if h is not affine, / is strietly inereasing every¬ 
where. To see the latter, suppose h! is not eonstant and fix any 6 G M. Then, for suffieiently small —M, 
we must have h''{y + b)dy + E-M<y<oo > 0 with A/i'(f/) := h'{y+) - h'{y-). 

Then, beeause h' is nondeereasing. 


I\b) > 


if{q) - <l)(g) 

</(7) 


h'\y + b)e-^^^'^ydy+ ^ A/i'(|/+ 


0<y<oo 


, nO poo _ poo . 

+ 5(^j h'\y + b) j e-^^^>Q^y\z-y)dzdy+ ^ Ah'{y + b)l - y)dzj 


-M<y<0 


> 0 . 


Therefore, for the affine ease with / = 0 (whieh holds if and only if h'{b) = q(3 uniformly in 6 G M by 
Lemma 3.2 below), we set b* to be any value on M. Otherwise, b* (if —00 < b* < 00 ) beeomes the unique 
root of /(/) = 0 . 

Example 3.1. For the case h{y) := ay'^, y G M., for some a > 0 , 6* = I3q/{2a) + E(-XgJ - ^{q) \ 

For the eases b* = 00 and b* = — 00 , we shall set our eandidate value funetions to be Vooix) := 
limsupb.|.oon 6 (x) and V-cx^ix) := liminf{, 4^.00 ^^(x), respeetively. By (2.8), if we set := inf{f > 0 : 
Xt > /jandfi:^ := inf{t > 0 : T) < /}, theestimation( 2 . 8 )immediately givesEa,[J°[?e“^*|/i(/ 4 ^)|df] 

'^b 
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II 


Oand 


b].— OQ 


> 0. Hence, we have probabilistic representations: 


(3.28) 


Fqo (t:) Ej 




and V-r^ix) = 




/36 

3-) 

Q 


which are the NPVs corresponding to the strategies i°° = 0 and i °° = By the resolvent measures for 
X and Y, we can write (3.28) succinctly using the scale function: by Corollary 8.9 of [14], 


(3.29) 


Vc^(X = 




Jo 

1 


h{x + z)e + / h{x + z)u^'^\—z)dz, 

J —CO 

[ h{x + z)e~'^^‘^^^dz + [ h{x + z)u^‘^\—z)dz + —, 

Jo J —OO 0. 


where := and := — W^'^^(t<;), for all tp > 0. 

The criteria for b* = oo and b* = —oo can be obtained concisely below. Let h'{oo) := limy_^oo h'{y) G 
(—cxo, cxo] and h'{—oo) := linij^^.oo h'{y) G [—cxo, cxo) (which exist by the convexity of h). 


Lemma 3.2. We have 


6<!>{q) fh'{oo) \ J T/ \ 5^{q) h'{-oc) \ 

ip[q) \ q J ip{q) \ q J 

Hence, I{—oo) = I{oo) = 0 if and only if h'{b) = q(3 for all b E M. Otherwise, b* = oo if and only if 
h'{oo) < q(3 and b* = —oo if and only ifh\—oo) > q(3. 


Proof By monotone convergence applied to the first equality of (3.23), 


I{oo) = h'{oo) 


ip{q) - $(q) 


^0 POO 


+ (5 


.-‘P(^p0(^)(z-y)dzdy 


' —OO J 0 


-P6 


Hq) 

piq)' 


9?(g)2 

Hence, when h'{oo) = oo, /(oo) = oo. 

Suppose h'{oo) < oo. By Fubini’s theorem, integration by parts and Remark 3.2 (1), we have 


^0 POO 

/ / ' 

'-oo Jo 

m f 

q Jo 

^ ^(q) 
q 


^-^^<ipQ(i){z-y)dzdy = - y)dydz 


3 -v(9 )^P|_X^ Z>z}dz = 


^{q) 

q 


r 1 _ p-v{q)^ 

' Gdz} -—— 

(0,oo) ^ ^{q) 


J (0,oo) 
q(p{q) 


[-X, edz} / 


1 - 


[0,oo) 


Gdz} 


^(g) n - Fp^(?)^eo = 

qpiqy ’ qp{q) ^{q?5 ’ 


where the last equality holds by (3.9). Substituting this, we have the expression for J(oo). The proof for 
J(—oo) is similar. □ 


The following example is a direct consequence of Lemma 3.2. 


Example 3.2. For the case h{y) := ay + y, y E M., for some a,r] E M., we have b* = —oo when a/q> (3, 
b* = oo when a/q > [3, and any value on [— oo, oo] when a/q = (3. 
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Using the obtained b* G [— 00 , 00 ], our eandidate optimal strategy is with its expeeted NPV given by 
Vh*. For the ease h* G (—cxo, cx)), by our ehoiee that ui,* (x) = 0, we have 


(3.30) 


Ufc* (x) = 


)dt 


x G M, 


'-Jo 

where the differentiability at h* holds beeause the right hand side is eontinuous dXx = h*. 

In view of (3.29), the relation (3.30) also holds for b* = —00 and b* = 00 . To see how the derivative 
ean go into the integral in (3.29), by the eonvexity of h, we ean ehoose a suffieiently small —m sueh 
that h{x + z) is of the same sign for all z < —m — x. Henee, the eonvexity of h implies by monotone 
eonvergence that 



\h{x + z)-h{x + e + z)\ ^^(^^^ 


e-s-O 
-)■ 


\h'{x + z)\u^'^\—z)dz] 


the same ean be said when is replaeed with 


4. Verification of optimafity 

We shall now show the optimality of the strategy and the assoeiated expeeted NPV of total eosts 
Vb*. Toward this end, we shall first obtain a suffieient eondition of optimality (Lemma 4.1 below) and then 
show that Vb* satisfies it. 

4.1. Verification lemma. For the given speetrally negative Levy process X, we call a function / suffi¬ 
ciently smooth if / is continuously differentiable (resp. twice continuously differentiable) on M when X 
has paths of bounded (resp. unbounded) variation. We let F be the operator acting on sufficiently smooth 
functions /, defined by 

2 r 

Tf{x) := yf{x) + ^fix) + / [f{x + z)- f{x) - /(x)zl{_i<^<o}]i^(dz), 

J (—00,0) 

and let (F — q)f{x) := Tf{x) — qf{x). If, for some i G Fl^, is sufficiently smooth, then Tve is 
well-defined. In addition, it is finite by (2.8) and Assumptions 2.1 and 2.2. 

Lemma 4.1 (Verification lemma). Suppose a strategy f G Fl^ is such that is sufficiently smooth on M 
and satisfies 

I (r - q)vfix) + h{x) > 0 ifv\{x) < 

I (F - q)vfix) - 5{v'^{x) - /3) + h{x) > 0 ifv'^^{x) > fi. 

Then I is an optimal strategy and v{x) = vfix) for all a; G M. 

Proof. We first note that (4.1) is equivalent to the condition 

(4.2) ^inf ^{(F - qffifix) - r(n'-(x) - fi)} + h{x) > 0. 

Hence, we assume (4.2). For brevity, we write w := throughout the proof. 
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Fix any admissible strategy t G and let (T„)„gN be a sequence of stopping times defined by T„ : = 
inf{t > 0 : |f//| > n}. Since = X—L is a semi-martingale and te is sufficiently smooth by assumption, 
we can use the change of variables/Ito’s formula (cf. [19], Theorems 11.31 and 11.32) to the stopped process 
t > 0} and deduce under P^. that 

rtATn rtATn 

= - / e-^^qwiUl)ds + / e-''V(f/f_)d(X, - L,) 

-'0 Jo 

2 rtATn 

+ — e-''W'(f/f_)ds + ^ + AX,) - AX,], 

0<S<tATn 

where AC, := C “ C- ^^d Aw((s) := w((s) — tF(C-) for any right continuous process (. Rewriting the 
above equation leads to 

rtATn rtATn 

- ifW = / e-»*(r - g)w(C/,'_)ds - / e-«-w'{U‘,_)AL, + 

Jo Jo 

with 

ae-'^W(f/f_)d5, + lim / [ e-'^^w\U^^_)y{N{ds x dy) - v(d|/)ds) 

'^'1° J{0,u] J{-l-e) 

+ [ [ + y)- w{Ul_) - w\Ul_)yl{y^(-i^Q)}){N{ds x dy) - z/(di/)ds), m > 0 , 

J (0,m] J (—cxDjO) 

where {i?,; s > 0} is a standard Brownian motion and X is a Poisson random measure in the measure 
space ([0,cxo) x (—cx), 0), S[0, cxo) x ;B(—cxo, 0), ds x z/(da;)). By the compensation formula (cf. [14], 
Corollary 4.6), t > 0} is a zero-mean martingale. 

Now we write 

rtATn 

w{x) = - / e-"* [(r - q)w{Ul_) - e,{w'{u‘,_) - P)] ds 

Jo 

rtATn 

+ P / e-^^^ds + - MtATn • 

Jo 

By (4.2), after taking expectation, w{x) < e“'?®(/i((7f) -|- /3£,)ds] + 

By (2.8), linij^nioo e“'?^/i(f/f)ds] = E 2 ;[Jg°° e“'^^/i((7f)ds] by dominated convergence. On the 

other hand, because i is nonnegative, we also have limj ,ntoo^4C^^ e-^'/3isds] = E^i/^e-^^/^Ads] via 
monotone convergence. Finally, by 0 < £, < 5, (2.6), (2.7) and the strong Markov property of X and Y, 

^ e“'^^(|/i(X,)| + |/i(X,)| + |/i|l{L>-oo} + /9<5)ds]. Therefore, Fatou’s lemma 

and (2.8) give limsup^^^^^ < 0. 

Hence, w{x) < Vi{x). Since £ G was chosen arbitrarily, we have that w{x) < v{x). On the other 
hand, we also have the reverse inequality, w{x) > v{x), because w is attained by an admissible strategy i. 
This completes the proof. 

□ 
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4.2. Optimality of Vb*. In view of Lemmas 4.1, it suffiees to show that the funetion Vb* is suffieiently 
smooth and satisfies (4.1). The former ean be verified by using the form of as in (3.30) and eonsidering 
its derivative; we defer the proof to Appendix B. 


Lemma 4.2. The function Vb* is sufficiently smooth. 

We shall now verify the inequalities (4.1). Toward this end, we restate it in the following lemma. 


Lemma 4.3. The inequalities (4.1)/or = Vb* hold if and only if 


(4.3) 


v'b,ix)>P ifx>b*, 
x'b*(x) < (3 ifx < b*. 


Proof, (i) We shall first prove that the following equalities hold for Vb*: 

(r — q)vb* {x) + h{x) =0 for X <b*, 
(r — q)vb* (x) — 5{v'b* {x) — /9) + h{x) = 0 for X > b*. 


(4.4) 


Here we give a proof of the seeond equality of (4.4); for the first equality a similar (and even simpler) 
argument holds. We also foeus on the ease b* E {—oo,oo); the proof for the eases b* = —oo and b* = oo 
ean be obtained by modifying the hitting time k defined below. 

Fix X > b* and reeall (2.4) for the definition of Y. Define the stoehastie proeess 

pt/\K 

Mb ;= + /3<5)ds, t > 0, 

Jo 

where n := inf{f >0:1)^ [b*, A^]} with N > x. 

In order to show the seeond equality of (4.4), it suffices to show that M is a P^;-martingale. Indeed, 
it holds by the martingale property and Ito’s formula (thanks to Lemma 4.2, which guarantees that Vb* is 
sufficiently smooth) noticing that the generator of Y is given by Ty/ := T/ — 5f'. 

For X > b* and t > 0, the strong Markov property of Y gives 




e-'^^Vb4Yn) + / e-'?*(h(y,) -f f35)ds 




= h 


{Kk} ) ® 


e-^^Vb4y.)+ e-‘^^{h{Y,)+P6)ds 


+ / e '^^{h{Ys) -f /3(5)ds > + 


Because U^* = Yg for all 0 < s < k, P 3 ;-a.s., 


d{t<K}Mt =l{t<K} < e 


e-^"nfe*(r,)+ / e-'^^{h{Yg)+(35)ds 


+ / e-'^^{h{Yg)+(36)ds\. 


Putting the pieces together, we conclude that Mt = Ej; [e~'^'^Vb* (Yfj + e~'^^(h(Yg) -|- /35)ds| ^ > 0, 

which is a martingale. 

(ii) Using the identities (4.4), we shall now complete the proof of this lemma. 

Suppose first that (4.3) holds. If u (x) < fi, then (4.3) implies x < b* and so, by (4.4), (F — q)vb* (x) + 
h(x) = 0. If f(,.(x) > (3, then (4.3) implies x > b* and so, by (4.4), (T— q)vb*{x)— 6{Vb,{x) — (3)+ h{x) = 
0. If fb*(a;) = [3, then we have, by (4.4), (F — q)vb*{x) + h{x) = 0. Hence (4.1) holds with = Vb*. 
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To conclude the proof, suppose now that (4.1) holds with vi = Vb*. For the case x < b*, suppose for 
eontradietion that > {3. Then (4.1) and (4.4) give —6{v'fj,{x) — (3) > 0 whieh implies vl,{x) < (3 

and henee forms a eontradietion. Therefore we deduee v'^, (x) < [3. Similarly, for the ease x > b*, suppose 
u^*(a;) < (3. Then (4.1) and (4.4) give 6{v'fj,{x) — (3) > 0 whieh implies ^^*( 0 ;) > (3 and henee forms a 
eontradietion. Therefore we deduee Ffe* {x) > 13. 

□ 


Lemma 4.4. The function Vb* is convex. 


Proof. By (3.30), it is suffieient to show that the mapping x i-A is nondeereasing. 

(i) Suppose that — 00 < 6* < 00 . Notice that the law of Uj’* under Pa- is the same as under P where 
U\ is the unique solution to the SDE: '^ = x + Xt — 5 f* By the convexity of h its 

derivative h' is nondeereasing, and henee it is suffieient to show that, for any fixed y > X, Uf’y > f/f’" 
P-a.s. for all f > 0. 

First, let us suppose that X is of bounded variation. We define a sequenee of inereasing stopping times 
0 < Ti < T 2 < ■ ■ ■ representing the times at which either or crosses the level b*. For example, 
ify>x>b*, then ti = inf{f > 0 : ’^ < b*} and T 2 = inf{f > ti : Uj’ < b* or > b*}. 

By induetion, we shall show that for any t G r„] for eaeh n > 0 with tq = r_i ;= 0 

by eonvention. The base case (n = 0) is already clear, because 'V = y>x = 

Now suppose that the hypothesis holds for some n > 0; this implies . Then, there are 

three possible scenarios (1) Ufj > > b*, (2) U^y >b*> Uy \and (3) b* > yy > yy. 


For the ease (1), it is neeessary that = yy + {Xt — xy — 6{t — Tn) > yy + (Xt — xy — 

b(t- Tn) = y for Tn<t < Tn +1 = inf{f > : y < b*}. 

For the ease (2), we have r„+i := Tn+i A r„+i, where r„+i := inf{t > Tn : y < b*} and 

Zn+i := inf{t > Tn : > b*}. ^ 

On the event {r„+i < r„+i}, > b* > f/ffor all f G [r^, Tn+i) and yy^ = yy^_+AXr„^^ > 

yy_ + AXn^^, = yy. On the event {r.+i > J, f/f’^ >b*> f/f’" for all t G [r„, r„+i]. 

For the ease (3), we have f/f= yy + {Xt - Xy > yy + {Xt - Xy = f/ffor Tn < t < 


Tn+l = inf{t > Tn : yy > b*]. 

Henee for all oases y > y for all t G (r„, r„+i]. By mathematieal induetion, the inequality holds 
for all f > 0 . 

For the case of unbounded variation, recall, as in Lemma 12 of [15], that there is a sequence of refraoted 
Levy proeesses of bounded variation that eonverges a.s. to the desired refraeted Levy proeess. Henee, we 
oan obtain the same inequality y > y P-a.s. by taking the limit. 

(ii) Suppose b* = 00 or b* = — 00 . In view of (3.28), beeause Voo{x) = E[J“ e~'^^h{Xt + x)dt], the 
eonvexity of h gives v'yx) = E[jy e~'^^h'{Xt + x)dt] = e~'^^h'{Xt)dt]. Similarly, v'_yx) = 

E^ljy e~'^^h'{Yt)dt]. The eonvexity of h now shows that x 1 —)■ Ej; [jy e~'^^h'{y*)dt\ is nondeereasing. 


□ 


Proposition 4.1. The function Vb* satisfies {A.\). 
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Proof. In view of Lemma 4.3, we show that Vh* satisfies (4.3). 

(i) Suppose —cxo < b* < oo. By setting x = b* in (3.30), we have 

piq) - 


v[,{b*)= / h'{y + b*)- 




+ 


h\y + h 


5^{q) 

j- V7(g) -$(g) 

n $(g) 




= /( 6 *)^^ + /3 = / 3 . 


(5$(g) 

This together with Lemma 4.4 shows (4.3). 

(ii) Suppose 6* = oo or 6* = — oo. For the case b* = oo, by the convexity of h, monotone convergence 
gives lim 3 ;^oo= lim^too e“^*h'(JA't + x)dt] = h'{oo)/q < /3, where the last inequality holds 
by Lemma 3.2. Similarly, for the case b* = —oo, lim 2 ; 4 ,_oo v'_^{x) = lima; 4 ,_oo E[/q°° + x)dt] = 

h'{—oo)/q > {3. These bounds together with Lemma 4.4 show (4.3). 

□ 


By Lemmas 4.2 and 4.3 and Proposition 4.1, we now have the main result of the paper. 

Theorem 4.1. The strategy is optimal and the value function is given by v{x) = Vb* (x) for all a: G M. 

We conclude this section with the remark on how the convexity of h is important in deriving our results. 

(1) Note that the function I as in (3.23) is monotone due to the convexity of h. Otherwise, there 
may be multiple b* such that I{b*) vanishes; the form of the optimal strategy will be much more 
complicated. 

(2) Thanks to (3.30) and the monotonicity of h!, is also monotone. This step is necessary in 
showing the sufficient condition for optimality (4.3). 

(3) The convexity of h also helps in technical details; indeed by the monotonicity of h, we can ex¬ 
change limits over integrals in various parts. 

(4) The uniqueness of the optimal strategy t'* is not guaranteed in the complete set of admissible 
strategies II^, but it holds when we restrict to the set of refracted strategies when h is not affine. 


5. Convergence to Reelection Strategies 


Recall from Remark 2.1 that 


(5.1) 


vix] 5) := inf E^. 

IdUs 


+ =v{x-5,-P) + 


/3(5 


-JO - q 

where n(x; 5, —/3) is the value function (2.3) obtained in the previous sections with Xt replaced with 
:= Yt + 5t and f3 with —/3. In this section, we fix the process Y (with its Laplace exponent V^y(-)> 
p{q) := sup{A > 0 : (A) = q] and scale function W(-)) and then study the asymptotics as 5 ^ oo- 

When the absolutely continuous assumption on L is relaxed, the problem defined in (5.1) becomes 
a classical singular control problem. Let Iloo be the set of admissible strategies consisting of all right- 
continuous, nondecreasing and adapted processes L with Lq- = 0. It is known as in Yamazaki [21] under 
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Assumption 5.1 defined below that, for all x G M, 


vix: oo : 


: = inf 




'[ 0 ,cxd) 


e + Lt)(it + /3dLt) 


(5.2) 


n i \\ L V'y(0+) 

= —p Z [x — b (cxo)) + 


Q 


' 6*(oo) 


— y)h{y)(\.y 


+ Z^^\x-h*{oo)) 




q Jo 


M,)yh(y + b*{oo))dy+^), 


where := 1 + q W^'^\y)dy and z'''^\z) := Z^‘^^y)dy for all 2 ; G M. The infimum is attained 

by the reflected Levy process Yt + with 

pb ( 00 ) [[b*{oo)) — Yt') VO, t > 0. 

0<t'<t 


The lower boundary 6* ( 00 ) is defined as the unique root of Joo(6) = 0 where 


(5.3) 


loo(b) : = 



h'(y + + /3 




b eR. 


Our objeetive in this seetion is to show the eonvergenees of b*(b) to b*(oo) and v(x;b) to ^(a:; 00 ) as 
5 t 

Throughout, we assume that F is a speetrally negative Levy proeess and satisfies Assumption 2.1 (2), 
whieh means that there exists 0 > 0 sueh that exp flyi—O) = E[exp(—0Yi)] < 00 . In addition, we further 
assume the eondition postulated in Yamazaki [21]: 


Assumption 5.1. (1) For some a G M, the function h{x) + fJqx is decreasing on (—cxo, a) and increasing 
on (a, cxo). (2) There exist a cq > 0 and an xq > a such that h'{x) + (Jq> cofor a.e. x > xq. 


We use the explieitly obtained expression of (5.1), for eaeh 5 > 0, with 

/3 = -/), 

and take the limit as 5 t oo- order to emphasize the dependenee on 5, let us denote, for all 5 > 0, 
flfls) := flyis) + 5s as the Laplaee exponent of the speetrally negative Levy proeess The seale 
funetions W'f'* and 0^^^ are defined in an obvious way. Also, let $ 5 (g) := sup{A > 0 : = q} and 


(5.4) 

h{h) := ~ r h\y+b)e-^^‘^^ydy+5 

Tiq) Jo 


h'{y+b) [ e-'^^y^^e^-‘\z-y)dzdy-/3^4^ 


Tiq) J 


Assumption 5.1 guarantees that h'{—oo)— jJq < 0 and h'{oo)—(3q > 0. Henee, by Lemma 3.2, Is{oo) > 0 
and Is{—oo) < 0, implying that there always exists a root b*{6) G (—C) 0 , 00 ). As in the diseussion given 
in Seetion 3.4, Assumption 5.1 guarantees that the root b*{5) is unique. We also let the 

eorresponding optimally eontrolled proeess. In other words, it is a refraeted Levy proeess defined as the 
solution to the SDE: duj;^^’^ f > 0. 
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Note that we have $ 5 (g) 0; beeause q = i/jyi^siq)) + and xfjy is eontinuous on [0, oo) and 

vanishes at zero, 


(5.5) 


<5<D,(g) ^ q. 


By (8.20) of [14], for any 9 > ip{q) (whieh is uniformly larger than ^s{q)), /poo)® ^^5W'f\(lx) = 
sd/iipyie) + 5e-q) 

(5.6) 


A 1. Henee continuity theorem gives that 


(x) = [0, x] ^ 1, X > 0. 

5.1. Convergence of optimal thresholds. We first show the convergence of b*{5) to b*{oo). Toward this 
end, we first show the following. 


(Ifoo 


Lemma 5.1. For all 6 G M, h{b) - Ioo{b). 

Proof. It is clear that /g°° h’{y + dy //° h’{y + Hence, it is left to 

show that //^ \ h'{y + b) \ //° — y)dzdy 0. 

Recalling Assumption 2.1 (2), we choose 0 < e < 0. For any y < 0, 

poo poo poo 

/ e-^^'^>5e^^\z - y)dz = 


'-y 


' -y 


which is bounded from above by exp(e|/) times a constant because (with and Y_ the running infimum 
processes of X^^'> and Y, respectively) 

[ e^^6el^\w)dw < 


-E[e 


=ireJi 


'-y 


which is bounded in 5 by (3.11) (with X replaced with Y) and (5.5). This, together with //^ \ h'{y + 
bfe^'^dy < cxo by Assumption 2.2, allows us to apply the dominated convergence theorem and 


lim 

5—)-oo 


/•exj p\j /‘Lxj 

\h\y + b)\ e-^^'^^^6ef{z-y)dzdy= \h'{y + b)\\im {z)dzdy. 


5—>-oo 


-y 


Again, by Remark 3.2 (1) and (5.5), for any 9 > 0, 


(■5), 


E[e^Aeq ] ^ 


q S{9-^s{q)) 5too, 


> 1 . 


S^s{q) i>Y{9) + 69 - q 
Hence Xj ^^ converges to zero in distribution as <5 t oo- Using this and (5.5), for all y < 0, 

f < f 6e^^\z)dz = ^ 0 . 


'-y 


' -y 


This shows the claim. 


□ 


By Lemma 5.1, we now show the convergence of b*{5). 


Proposition 5.1. We have b*{5) 6*(oo). 
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Proof. Let 5 > 0 be large enough such that ^sio) < 1/2 for all 6 > 6. 
By differentiating (5.3), 



b eR. 


In addition, by (5.4) and because h' is nondecreasing, we have a lower bound: 
(5.8) heR,5>5_. 


Notice that Ioo{b) is strictly increasing on (—oo,6) for some b E (— 00 , 00 ]. When 6 < 00 , we must 
have looib) > 0; otherwise Assumption 5.1 (2) is violated. This means that we can choose some b G 
(6*(oo),6) such that Ig{b) > 0. In addition, for all b < b we can obtain from (5.7) a lower bound 
^'ooib) > r^(b). Hence we can choose a sufficiently small e > 0 and a > 0 such that I'^{b) > a 

on {b*{oo) — e, b*{oo) + e); this together with (5.8) gives 


(5.9) 


inf 

S>5, b^{b* {oG)—£,b* (oo)+£) 


m > f > 0. 


Fix any 0 < e < A By Lemma 5.1, we can choose a sufficiently large 6 > 6 such that |Joo(&*(oo)) — 
Is{b*{oo))\ = |/ 5 ( 6 *(oo))| < ea/2 for any 6 > 6. Then, by (5.9), we must have for any 6 > 6 that 
Is{b*{oo) — e) < Is{b*{oo)) — ea/2 < 0; similarly, Is{b*{oo) + e) > Is{b*{oo)) + ea/2 > 0. This 
means, together with the monotonicity of Is, that b*{oo) — e < b*(S) < b*(oo) + £. Because e was chosen 
arbitrarily, this completes the proof. 

□ 


5.2. Convergence of value functions. With the help of Proposition 5.1, we show the convergence of 
v(x; S) to v(x; 00 ). More precisely, we show the following theorem. 

Theorem 5.1. Uniformly in x in compacts, v{x] 6) v{x] 00 ). 

In order to show this theorem, it is sufficient to show the pointwise convergence. Indeed, by the defini¬ 
tion of v{x-, (5) in (5.1) as the infimum of the NPV and because the set II^ is monotonically nondecreasing 
in 5, v{x; 5) is nonincreasing in 6. Hence, pointwise convergence implies the uniform convergence on 
compacts by Dini’s theorem. 

We shall first show for x < 6* ( 00 ) and then extend the result for x > b*{oo). 

Lemma 5.2. For any x < b*{oo), we have v{x] 5) v{x] 00 ). 

Proof. By Proposition 5.1, we can choose sufficiently large 5 such that x < b*{6) for all 5 > A Let 
h{—oo) := lim 64 ,_oo h{b) E [— 00 , 00 ], which exists by the convexity of h. 

Fix any e > 0. 

(i) Suppose h{—oo) = 00 . By the convexity of h (together with Assumption 2.2) and because (3.11) 
holds when X is replaced with Y, we can choose a sufficiently small —M such that h is decreasing and 
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positive on (—oo, — M] and 








= E, 




'-JO 


< e. 


Because Y_^ < a.s. for any S > S, 


(5.10) 


sup E^ 




< e. 


(ii) Suppose h{—oo) G (—cxo, cxo), then h is bounded on the negative half line and we can choose —M 
such that (5.10) holds. 

(hi) Suppose h{—oo) = —oo. Then, we can choose —M such that h is negative on (—cxo, —M] and 
hence *''^^)l.j^(d).6*(d)^_^,dt] < 0. In view of these, it is sufficient to show for all 

-M e M 


v{x; 5, —M) := v{x] 5) — E^ 


Jo 






> v[x; OO). 


Indeed, this implies for the cases (i) and (ii) that lim^^co |u(x; 5) — i;(a:; oo) | < £ for any e > 0. For the case 
(hi), it implies limsup^^oo v{x; 6) < v{x] oo) (and liminf^ioo v{x] 6) > v{x] oo) holds because C IIoo 
in view of (5.1)). 

Fix S > S. We can write 

v(x; S, -M) = r f ^ 5*(5))e-v^('?)yd|/ + ^ 

o^s{q) LJo 


+ 


'-M-b*{5) 


h{y + b*{6)) 


Mq)(x-b*(s)) ^iq) - ^s{q) 

Mq) 


(^(g)j q 
-y)dz- W^'^\x - b*{6) - y) 


dy. 


Regarding the first line on the right hand side, it equals 


(( 5 )) 


y?(g) - ^sjq) 
S^s{q) 


h{y + b*{5))e-^^‘i^ydy + P 


^s{q) q qj{q) 


^(1 _ Q^s{q){x-b*{5))^ 


(5too 


> v{x] oo) = (3 \x — b*{oo) + 


V^y( 0 +) \ ^ ^ 
q J q 

To see how the convergence holds, note, by (5.5), 


-^iq)yh{y + b*ioo))dy-^. 

^{q) 


(5.11) 

In addition. 


^.i.,(,)(x-b*(S)) ^iq) - ^s{q) stc^ ifjq) 
d^q) q • 


1 6 _ q- S^s{q) _ -ipYi^siq)) 5too^ '0y(O+) 


^siq) q ^s{q)q ^s{q)q 
and the Taylor expansion and the convergence of b*{6) gives 

5(1 _ e-*-d(q)(x-fe‘(5))) ^ 5 [(^ _ b*{6))^s{q) + o(<f>5(g))] ^ q{x - b*{oo)). 
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On the other hand, —M — b*{5) eonverges and hence is bounded in 5 from below by, say —N < 0; 
hence it is now left to show that 


(5.12) 

^<I>5 (q) {x-b* (S)) 9 ?(g) — ^sjq) 

Mq) 

Here, sup_^<y<os> 5 \h'{y + < oo and 


\h{y+b*m 


’-N 


,-^{q)z^^{qy^z-y)dz-W^'^\x-b*{6)-y) dy ^ 0 . 


fO , 


0 < 


(5.13) 


l-N 


<i>s{q){x-b*{S)) ^i^) ^s{(l) f (^ 2 : — y)dz 


^siq) 

- Mq)(x-b*(5)) ^iq) - ^s{q) ^ 


dy 


= e 


6^s{q) 

Here we recall the convergence (5.11). In addition, 

poo 

6 / + N)- W^'^\z))dz 

Jo 

r-N 


+ N)- W^‘^\z))dz. 


(5.14) = I e-‘^^'^'>^W^'^\z)dz- I e-^^'^'>^W^‘^\z)dz^ - 6 j e-^^^^^W^'^\z)dz 

pN 

- / e-^^'^'>^6W^'^\z)dz ^ 0 . 


— U'Piq)N _ 


1 )- 


5 


JjY{.v{.q)) + bip{q) -q 

To see how the convergence holds, 5 supg< 2 <Ar W^'^\z) = 5Wg'^\N), which is bounded in 5 by (5.6). 
Hence dominated convergence and (5.6) give lim 5 .|.oo e~'^^‘^^^dz = (1 — 
e-‘^('^)^)/v9(g). On the other hand, (x — b* (5) — I/) —)■ 0 uniformly on 1 / G [—N, 0] by the convergence 
of b*{6) and (5.6). Hence (5.12) holds, as desired. 

□ 

In order to extend the result to the case x > 6*( 00 ), we shall show that the derivative converges. By 
(3.30), 


(5.15) 


v\x-, 5) = E J / e-^*/i'((/f 


X G 


On the other hand, by [21], the derivative of v{x; 00 ) simplifies and 


v'{x]oo)=j3— / {h'{y) — f3q)W^^\x — y)dy^ a; G 


' b* (00) 


(5too 


Lemma 5.3. For every x G M, we have vjx; 5) —^—)■ vjx; 00 ) for all x 7 ^ b*{oo). 

Proof. Fix any 5 > 0. We first show, for any c > 0, that we can choose sufficiently small —M G M such 
that 


(5.16) 


sup E 3 ; 

S>5 Ljo 




><-My 


< e. 


(i) Suppose h'{—oo) < 0. By the convexity of h and (3.10), we can choose a sufficiently small —M such 
that, on (— 00 , —M], h! is negative (and hence \h'\ is decreasing) and Ej.[/g°° e“'^*|h'(F))|l{Y(<_M}df] < e. 
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By r < for all 5 > 5, we have (5.16). (ii) Suppose h'{—oo) > 0, then the eonvexity of h implies 

that h' is uniformly bounded by a eonstant on the negative half line. Again, by F < (5.16) holds 

for suffieiently small —M. 

For any e > 0, with —M sueh that (5.16) holds, the relation (5.15) implies 


sup 

< 5><5 


v'{x] 5) - 












< e. 


Henee, for the proof of this lemma, it is suffieient to show that for all small —Mg 


v'{x; 6, -M) := E^ / e-'^*/i'((7f 


'-./o 




df 


(5too 


> V [X] oo). 


By (3.13) and (5.15), we can write v'{x] 6, —M) = v^^^'{x] 6, —M) + 5, —M) where 


- y)dz - W^'^\x - b*{5) - y)] dy, 


6, -M) := j h\y + h*dy 

J-M-b*{S) L $<5(0') Jo 

v^^^'{x;6,-M) := ^ h\y + b*{5))^^e-^^'^^yMs{x;b*{5))-W^'^\x + b*{6)-y)^dy 

/ O . roo 

h'{y + b*{S))\Ms{x;b*{6)) / - y)dz 

-M-b*{ 5 ) Jo 


'b*{S) 


W^'^\x - z)W^'^^'{z -y- 6*((5))d^}d|/, 


where Ms{x-, b*{6)) := {(p{q) — $ 5 (g))e ^dQ)f>*{S) Q^s{q)z'^{q )_ z)dz. 

Regarding v'^^'>'{x] 5, —M), by the choice of b*{oo) such that (5.3) vanishes. 


.<i>,(q)(x-b*(S)) ^iQ) - ^s{q) 

6Mq) 


h'iy + fe*(( 5 ))e-^('?)M 2 / ^ ^ [ h'iy + 6 *(cx)))e-^('')M 2 / = 

q Jo 


The remaining term of v^^'>'{x] S, —M) vanishes in the limit. Indeed, we can choose sufficiently small 
—N < 0 such that —M — b*{6) > —N uniformly in 5 > 5. Clearly, sup_jv<jy<o, 5>5 \ h'{y + &*(<)))| < 00 , 
and we have shown that (5.13) vanishes in the limit by (5.14). Hence, 


(5.17) 


;(^)'(x;(5,-M) ^/3. 


It is now left to show the convergence v^‘^'>'{x] S, —M)l{j.>b*( 5 )} 00 , —M)l^a:>b*{oo)}- First, 

Proposition 5.1 gives l{3;>fe*(5)} l{3;>fe*(oo)} for a; 7^ b*{oo). Hence, it remains to show that < 5 , —M) 

converges to u‘^^)'(a:; 00 , —M) for x > b*{oo). 
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By Lemma 5.1, we ean set 5 large enough sueh that x > b*{5) for all 5 > 5. For any sueh <5 > 5, we 
have the eonvergenee 


h\y + - b*{6) - y)^dy 


f*x—b*{oo) 


W^‘^\z)dz-W^‘^\x-b* (oo) — 


j h'{y + b*{oo))^^e j 

= - [ ih'{y) - Pq)W^^\x - y)dy, 

J b* (cxo) 

where we used the ehoiee of b*{oc) that makes (5.3) zero. We shall now show that the remaining term 
vanishes in the limit and henee 


(5.18) 


L(2)'(a;; 5) ^ - / {^h\y) - Pq)W^^\x - y)dy. 

Jb* (oo) 


For any 6 > 6, Fubini’s theorem gives 


’-M-b*{5) 


{Ms{x]b*{6)) / e-^^‘^^^W^s’^^'{z-y)dz- W^'^\x - {z - b*(6) - y)dz 

'■ Jo Jb*{S) 

5|m(x; b*{S)) + M + b*{S)) - W^’^\z))dz 

pX 

W('')(x - z){W^'^\z + M) - W^'^\z - b*{6)))dz 


lb*{S) 


Here 5 + M + b*{6)) - W^'^\z))dz < 6 {W^'>\z + M + x) - wl'^\z))dz, 

whieh vanishes as 5 t cx) by (5.14). On the other hand, with 6* := inf 5>5 6*((5) > —oo, we have 
sup 6 *( 5 )<^<^ -z)= - b*{8)) < - k*) < oo, and 

f-x-b*{S) 


0<6 I iW^'^\z + M) - - fc*((5)))d^ = 5 

Jb*{S) 


dq) I 


{W^‘^\z + M + b*{5)) - Wl'^’{z))dz 


(9)/ 


(*x—b* 


< 5 I {W^'^\z + M + x)- Wl'^\z))dz, 


(9)/ 


whieh vanishes in the limit as 5 t oo by dominated eonvergenee and (5.6). Henee, (5.18) holds. 


□ 


Proof of Theorem 5.1. In order to eomplete the proof of Theorem 5.1, it is left to show that v{x-,6) 
v{x-, oo) for any x > b*{oo). 

Fix any x < b*{oo). We have F(x; 5) = F(x; 5) + v'(z; 5)dz for all x > t. By Lemma 5.2, we have 

F(x; (5) F(x; oo). In addition, for dXX x< z < x, by the eonvexity of v as in Lemma 4.4, v'{z] 5) < 
v'{x] 5), whieh is bounded in 5 by Lemma 5.3. Henee, Fatou’s lemma gives limsup^ioo ^ 

f^v'(z;oo)dz = F(x; oo) — F(x; oo). This shows limsupF(a:; 5) < F(a;;oo). Beeause F(x; 5) > 

v(x; oo), we have v(x; S) v{x-, oo) for any x > b*{oo), as desired. 

□ 
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6. Examples under phase-type Levy processes 


In this section, we focus on the case the processes X (and Y) have i.i.d. phase-type distributed jumps 
[1] of the form Xt — Xq = + aBt — 0 < f < oo, for some 7 G M and a > 0. Here 

B = {Bt;t > 0} is a standard Brownian motion, N = {Nt;t > 0} is a Poisson process with arrival 
rate k, and Z = n = 1, 2,...} is an i.i.d. sequence of phase-type-distributed random variables with 
representation {m,cx,T); see [1]. These processes are assumed mutually independent. The class of 
processes of this type is important because it can approximate any spectrally negative Levy process in law 
(see [ 1 ] and [ 10 ]). 

The Laplace exponent (2.1) of X is then a rational function and hence the scale function (and hence 
their resolvent measures (3.12) as well) admit analytical expressions; straightforward computation yields 
the optimal threshold level b* as well as the value function Vh* ^ 

As numerical illustrations, we focus on the quadratic case h{z) = z^, z ^ M. We numerically verify 
the optimality and then study the behavior with respect to 5. Throughout, let us fix g = 0.05. Lor the 
processes X and Y, we let a = 0.2 and k = 1 and, for the jump size distribution, we use the phase-type 
distribution given by m = 6 that approximates the Weibull random variable with parameter (2,1) (see 
[10] for the values of T and ck). 

We fix 7 = 5.5, (5 = 5 (and hence jy = 0.5). As it has been shown in Section 3.4, the function / is 
increasing. Hence, in general the bisection method can be applied to identify the value of b*. Using the 
obtained b*, we compute the value function Vb*{x) = v[P(x) -f (x)l{ 3 ;> 6 *}. 

In order to confirm the optimality of the refraction strategy, we plot on the top of Ligure 1, for the cases 
(3 = —5 and (3 = 5, the value function V},*{x) along with suboptimal NPV’s Vb{x) for b e {b* — 1, b* — 
0.5, b* -f 0.5, 5* + 1}. We observe that Vb* is indeed minimal (uniformly in x) and b* is the sole choice 
that makes the value function smooth and convex. In addition, the slope at b* of Vb*{x) equals (3. This is 
consistent with our observation that Vb* (x) satisfies (4.3). 

We then analyze the convergence as (51 cxo. As we have studied in Section 5, we fix the process Y with 
7 y = 0.5 and confirm the convergence v{x; (5) of (5.1) to v{x-, 00 ) of (5.2). Throughout, we fix /3 = 5 
(P = —5). At the bottom of Ligure 1, we first plot in the left panel the values of b*{6) with respect to 6 
along with the value of b*{oo). We see that b*{6) is decreasing in b and indeed converges to b*(oo) as b 
gets large. This is consistent with Proposition 5.1. In the right panel, we show the value functions v(x;b) 
as in (5.1) for b G {1, 2,..., 20,40, 60, 80,100} along with the function v(x; 00 ). As has been shown in 
Theorem 5.1, v(x; b) is decreasing in b and converges to v(x; 00 ). 


Appendix A. Proof of Lemma 3.1 


A.l. Proof of (3.25). (i) In order to show that the derivative can go into the integral, we shall show that 
for sufficiently small —M, 


(A.l) 



Analytical expressions of these can be found in a supplementary note given in https://sites.google.com/siteAyamazak. 
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Figure 1. (Top) Plots of Vb{x) for the cases /3 = —5 (left) and /3 = 5 (right). Each panel shows Vb* (x) 
(solid) in comparison to Vb(x) (dotted) for b € {b* — 1, 6* — 0.5, b* + 0.5, b* + 1}. The point {b* ,Vb* (b*)) 
is indicated by the square mark on the solid line. Up-pointing and down-pointing triangles show the points 
{b,Vb{b)) for b > b* and b < b*, respectively. (Bottom) Plots of convergence as <5 f oo. The left panel 
shows b*{6) for S running from 1 to 100. The value of b*{oo) is indicated by the dotted line. On the right 
panel, the functions v{x] <5) as in (5.1) are shown as dotted lines for (5 G {1, 2,..., 20,40,60, 80,100}. The 
solid line shows the function v(x; oo). The square mark indicates the point (b*(oo), v(b*(oo); oo)) while the 
up-pointing triangles show {b*{S), v{b*{S); (5)). 


We first see that, for e > 0, 

{x + e,y) = {x, y) — {x + e — y) + [x — y), y < b, 


and hence 

,( 1 ) 


rl"\x + e,y) -r^^\x,y) 




,( 1 ) 


(a;, y) _ + e - I/) - - y) 


e e e 

Here — is bounded in e > 0 on compacts and (x) = i/)di/is well-defined. 

Hence, for the proof of (A.l), it is sufficient to show that 

jx + e-y)- W^(^g) {x - y) 


y ^ H{y,x,e) : = 


h{y)e~ 


is bounded in e > 0 by a function that is integrable over (—cxo, —M). 
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First we observe that, by Remark 3.1 (3), (kF$(g)(x + e — ?/) — kF$(q)(a: — |/))/e < supo<g/<ekF4(q)((a:+ 
e' — y)+). Moreover, for any y < 0 and 0 < e' < e. Remark 3.1 (3) gives 

- y)+) ^ kF$(g ){x + e' - y) ^ + e - y) 

^4(g)((^-2/)+) “ W^{q){x-y) - W^^q){x-y) 

Fix any c, /c > 0. By (3.5) (and henee kF$(g)(a; + k)/Wq,i^q){x) 1), we ean ehoose suffieiently small 

— M < 0, sueh that for any y e (—cxo, —M], PF$(q)(a: + k — y)/Wq,i^q){x — y) < 1 + c. Colleeting these 
inequalities, we have a bound H{y;x,e) < (1 + — y)+)\h{y)\ for 0 < e < k and 

y < —M. Because (3.10) is finite, we have the claim. Hence (A.l) holds. 

(ii) By (i), the derivative can be interchanged over the integral. Hence, 


(A.2) 




^h{x)W^^\f)) + r f h{y + + 


/3(5 


+ / h{y + - ^{q)) / - y)dz - W^'i^'{x - b - y)\dy, 

J —oo ^ Jo 

where the first term appears due to the discontinuity of the scale function at zero as in Remark 3.1 (2) for 
the case of bounded variation. 

In order to apply the integration by parts, we first confirm that 


(A.3) h{y + b)\e'^^^^^^-'’\ip{q) - d>{q)) / - y)dz - W^‘^\x - b - y) 

L Jo 

Indeed, by (3.18), limy|_oo h{y + b)rl^\x, y + b) = 0. In addition, 

/•OO 

rl,^\x,y + b)- - d>{q)) / - y)dz - W^'^\x - b - y) 


yi-oo, 


> 0 . 


_ $(q)(x-b) y^(g) -^(g) r 

m Jo 

e < y} 

Using identity (3.11) together with Assumption 2.2, 
(A.4) 


- y)dz = <^}] 


n2L,<y}h{y + b)^^0. 


Hence, (A.3) holds. 

Using (3.2), (3.3) and (A.3), together with the following identity 


(A.5) 


i-^i<i>W^<i){z)dz = {'ip{ip{q)) - q)-^ = {^{q)b) 


-1 
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integration by parts yields 

Q-v^{<i)yh'(y + b)dy + P5 + h{b)W^'^\x - b) 

\ 

noo 

+ / h\y + b) -*^(?)) / - y)dz - W^'^\x - b - y) dy. 

J —oo ^ Jo 

Taking the differenee between this and h'{y)r\^\x, y)dy, we have the result. 

A.2. Proof of (3.26). In view of (3.17), let us deeompose vf\x) = vf'^\x) + 5vf''^\x) where 
yf’i)(a;):= [ (h{y + b)+/36)\e-‘^^‘>^yM{x-,b)-W^‘^\x-b- 


/ O ^ poo px 

h{y + b)^M{x]b) j - y)dz - I ^^^\x - z)W^^'^'{z - b - y)dz 


For the funetion vf'^\ integration by parts gives 


vf'^\x) = 


M{x]b) 


f'X—b 


w('^H2)d^|+ / h’ 


V>iq) 


fx—b—y 


W(«)(2)d;2}dl/. 


Differentiating this and noting that 

(A.6) M'{x; b) = ((^(g) - $(g))W('')(x -b) + $(g)M(a;; 6), 

we have 


(A.7) vf'^^\x) = 




{h(b) + 135) I {x-b) + M{x] b) 


+ 


+ t){ - ■K9))W<A^ - h) + ^ 


In order to deal with the funetion we shall first show the limit 
(A.8) 

where we define 

poo px 

N{y;x,b):=M{x;b) e-^^'i^^W^'^\z-y)dz- W^'i\x - z)W^'^\z - b - y)dz, 

Jo Jb 

so as to apply the integration by parts. 

First, beeause is well defined and finite by (3.19), limj^ 4 ,_oo h{y + b)r^^\x, y + b) = 0. In addition, 

poo px 

5-^/^\x,y+b)-d>{q)N{y;x,b) = M{x;b) / e-A'?)"0(‘')(2-r/)d^-/ W^'^\x-z)Q^'i\z-b-y)dz, 


lim h{y + b)N(y; x, b) = 0, 

j/T-00 
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whose absolute value is bounded by 

roo px 

M{x;b) / e-^^'^^^e^'^\z-y)dz+ / - z)e^'^\z - b - y)dz 

Jo Jb 

poo poo 

<M{x;b) / e-^('')"0(''H^-2/)dz + W^(,)(x-6) / - y)dz 

Jo Jo 

= <j,}] + W^(q)(x - <y}])^ 


< [M{x; b) + W^(g)(a; - < y} 




where W^(g)(z) := This together with (A.4) shows the limit (A.8). 

Now, integration by parts followed by Fubini’s theorem, together with the identities (A.5) and {x- 

z)W'^‘^\z — b)dz = W'''^^{z)dz), gives 




f*x—b 


W^'^\z)dz]+ W^'^\x-z)A{z;b)dz. 


Here we define 

A{z;b) := f h'{y + 


- y)dw - W^'^\z - y 


which is finite for all z e [b,x]. Indeed, we have &) = A{x-,b)+J^^h'{y+b){e~^^'^^^^~^'>W^'^'>{x- 

y — b) — W^'^^z — y — b))dy where the integral part is finite by similar arguments as in Remark 3.3; if it 
is infinity, it contradicts with that is finite. 

Hence, taking a derivative and by (A.6), 

vf'^^'ix) = - ^^SjpR + W^<i\x -b)- W^'^\x - 6)1 

6 I (p{q) J 


+ 


h'{y + b)B{y,x,b)dy 


^{^^(w^‘^\x-b)-M{x;b)^ -W^'^\x-b)^+ j h\y + b)B{y,x,b)dy 


with 

B{y,x,b) := 


- J - y)dw - W^‘^\x - y - b)^ 

/ X /* oo 

- y)dw - W^'^\z -y- b)^dz. 


By (3.24) and integration by parts, we can write 

{ifiq) - $(g))W('?)(x -b) + <I>(g)M(x; b) 


B{y,x,b) = 




- y)dz + W^^\-y) 


- W^‘^\x - b)W^‘^\-y) - / - z)W^‘^^\z -y- b)dz. 

Jb 

From this and (A.7), the expression J^^h'{y)rl^\x,y)dy — {vl^’^^'{x) + dvj^^’^^'(x)) gives the result 
after some calculations. 
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Appendix B. Proof of Lemma 4.2 


In view of the expression (3.30), is eontinuous for all a: G M. Henee, it remains to show that is 
eontinuously differentiable for the ease of unbounded variation. 

By (3.12) and (3.30), we can write where 


(5$(g) J 
+ f h\y + h*) 


h'{y + 


w^'^Hx) := 


' —00 
^00 




- y)dz - -b* -y) 


dy, 


h'{y + 6 *)|e-^('^)^M(a;; h*) - 


+ 6 h'{y + b*)\M{x;b*) - y)dz - W^‘^\x - z)W^'^^'{z - b* - y)dz 


Notice that (A.l) holds when h is replaced with h' (by Assumption 2.2). Hence, similarly to (A.2) 
(noting that = 0 for the case of unbounded variation). 


wd'>'(x) = 


ip{q) - $(g) 

5 


h'{y + 


+ / /i'(i/ + 6 *)p(^)(^-'’‘)((/?(g)-$(g)) / e-^^'^'^^W^'^^'{z-y)dz-W^'^'^'{x-b*-y)]dy. 

J—00 ^ Jo 

In addition, with the help of Fubini’s theorem and noticing that 1L*''^^(0) = 0, 


(2)'(a:) = h\y)b*) - W^'^^'{x - y)^dy 

/ X nO . roo 

W(^)'(a:-z) j /i'(i/ + &*)|((/?(g)-$(g))e-'*’(^)(^*-^) j - y)dw 

- W^'^^'iz -b* - 


The continuity of wdy and are clear for x 7 ^ b*. We have 

j h\y + fc*)|(</?(g) - <I>(g))e-'*’('^)(^*-^) ^ - y)dw - W^‘^^'{z - b* - 

= j h'{y + b*)i^{q)rd\z-b\y)-Q^‘i\z-b*- 

J —00 ^ 

By (3.7) and Remark 3.3, this is finite for 2 ; G [b*, x\. Hence, limj; 4 ^b* w^d'l^x) = 0. 
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